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Abstract. Using observed star formation rates at red- 
shifts up to z ~ 5, we calculate cosmic supernova rates 
for core collapse and Type la supernovae. Together with 
supernova statistics and detailed light curves, we estimate 
the number of supernovae, and their distribution in red- 
shift, that should be detectable in different filters with 
various instruments, including both existing and future 
telescopes, in particular the NGST. 

We find that the NGST should detect several tens of 
core collapse supernovae in a single frame. Most of these 
will be core collapse supernovae with 1 <,z < 2, but about 
one third will have z J> 2. Rates at z > 5 are highly uncer- 
tain. For ground based 8-10 m class telescopes we predict 
^0.1 supernova per square arcmin to Iab = 27, with 
about twice as many core collapse SNe as Type la's. The 
typical redshift will be z ~ 1, with an extended tail up 
to z ~ 2. Detectability of high redshift supernovae from 
ground is highly sensitive to the rest frame UV flux of the 
supernova, where line blanketing may decrease the rates 
severely in filters below 1 /jm. 

In addition to the standard 'Madau' star formation 
rate, we discuss alternative models with flat star forma- 
tion rate at high redshifts. Especially for supernovae at 
z > 2 the rates of these models differ considerably, when 
seen as a function of redshift. An advantage of using SNe 
to study the instantaneous star formation rate is that the 
SN rest frame optical to NIR is less affected by dust ex- 
tinction than the UV-light. However, if a large fraction 
of the star formation occurs in galaxies with a very large 
extinction the observed SN rate will be strongly affected. 
An additional advantage of using SNe is that these are not 
sensitive to selection effects caused by low surface bright- 
ness. 

Different aspects of the search strategy is discussed, 
and it is especially pointed out that unless the time inter- 
val between the observations spans at least 100 days for 
ground based searches, and one year for NGST, a large 
fraction of the Type IIP supernovae will be lost. Because 
of the time delay between the formation of the progeni- 
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tor star and the explosion, observations of z J> 1 Type la 
supernovae may distinguish different progenitor scenarios. 

A major problem is the determination of the redshift 
of these faint supernovae, and various alternatives are dis- 
cussed, including photometric redshifts. In practice a re- 
liable classification based on either spectroscopy or light 
curves requires the SNe to be ~ 2 magnitudes above the 
detection limit. The uncertainties in the estimates are 
discussed extensively. We also discuss how the estimated 
rates depend on cosmology. Finally, some comments on 
effects of metallicity are included. 



Key words: Cosmology: observations - Supernovae: gen- 
eral - Stars: formation - Nucleosynthesis 

1. Introduction 

The cosmic star formation rate (SFR) has now been esti- 
mated up to redshifts z ~ 5. By combining the evolution 
of the 2800 A luminosity density calculated from CFRS- 
galaxies (Lilly et al. 1996) at redshifts z < 1, and the 
1500 A luminosity density calculated from high redshift 
galaxies (z > 2) in the HDF found by the Lyman dropout 
techniques, Madau et al. (1996) argue that the SFR should 
peak at 1 < z < 2. Complementary to these, Connolly et 
al. (1997) have used HDF photometric measurements, to- 
gether with ground based near-IR photometry, to derive 
the 2800 A luminosity density at redshifts 0.5 <i z <i 2. 
These results are in good agreement in the region over- 
lapping with the CFRS, supporting the case of a peak 
in the SFR. This conclusion is, however, sensitive to the 
fact that the rest frame light may have been strongly at- 
tenuated by dust. Absorption by dust is especially severe 
for the UV-light, and since shorter rest frame wavelengths 
are sampled at higher z, this uncertainty increases with 
redshift. 

The evolution of the SFR is reflected in the cosmic su- 
pernova rate (SNR). It should therefore, in principle, be 
possible to use supernova (from now on SN) observations 
to distinguish between various star formation scenarios. 
Even more important, core collapse SNe, i.e. Types II 
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and Ib/c, provide a direct probe of the metallicity pro- 
duction with cosmic epoch. In reality these relations are 
non-trivial to establish. When it comes to core collapse 
SNe, the observational constraints at high redshifts are 
nearly non-existent. Even at low redshift the statistics 
are severely affected by selection effects. A main problem 
comes from the fact that core collapse SNe observationally 
show a large diversity, both in terms of luminosities and 
types, and with uncertain distributions. In addition, dust 
absorption, as well as background contamination, affect 
the statistics. Nevertheless, because of their importance 
for the nucleosynthesis, as well as galaxy formation, di- 
rect observations of the rate of core collapse SNe are of 
high interest. It is therefore hardly surprising that this 
is one of the main goals for the Next Generation Space 
Telescope (NGST) (Stockman 1997). 

For Type la SNe, the unknown time delay between 
formation and explosion of the progenitors unties the link 
to the SFR, making predictions more model dependent. 
Observations of the Type la rate at high redshift therefore 
provides a possibility to distinguish different progenitor 
scenarios. 

In this paper we present estimates for the expected 
number of observable SNe for the NGST, as well as for 
ground based instruments, and discuss various compli- 
cations entering the analysis. Previous studies include 
Madau et al. (1998a), Ruiz-Lapuente & Canal (1998), 
J0rgensen et al. (1998), Miralda-Escude & Rees (1997), 
Sadat et al. (1998), Yungelson & Livio (1997). With re- 
spect to most of these, our work differs in that we include 
information about the light curve, as well as spectral evo- 
lution, which allow us to predict the simultaneously ob- 
servable number of SNe. That this is important is obvious 
from the fact that a nearby SN seen at the tail of the 
light curve is indistinguishable from a more distant ob- 
ject at the peak. We divide the SNe into different types 
with maximum absolute magnitudes and spectral distribu- 
tions that varies with time and type. This also introduces 
a large dispersion in magnitudes at a given redshift. Ne- 
glect of these effects introduces a severe Malmquist bias. 
We calculate the counts for different broad band filters, 
and include information about the expected redshift dis- 
tribution of the detected SNe. Some preliminary results 
were given in Dahlen & Fransson (1998). 

Sect. 2 describes our model. Results are presented in 
Sect. 3. In Sect. 4 we discuss alternative star formation 
scenarios. In Sect. 5 we discuss how other cosmologies af- 
fect out results. The effects of gravitational lensing are 
discussed in Sect. 6. Problems concerning redshift deter- 
mination are discussed in Sect. 7. A general discussion 
follows in Sect. 8, and conclusions are given in Sect. 9. 
Throughout most of the paper we assume a flat cosmol- 
ogy with Hq = 50 km s" 1 Mpc 3 and Hm = 1, unless 
otherwise stated. 



2. The Model 

2.1. Filters and magnitude system 

Because of the strong UV deficiency in the spectrum of 
most SN types, filters bluer than R are of little interest 
for high-z studies. For R (A = 0.65 jum) and I (A = 0.8 
fim) we use Cousins filters, and in the near-IR J (A = 1.2 
fim) and K' (A = 2.1 /an) filters. For the M band we use 
a filter that we denote by M', which is centered on A = 
4.2 fim with A/AA = 3. For the magnitudes we use the 
AB-system where m A B = -2.51ogF„ - 48.6 (Oke & Gunn 
1983) {F v in ergs cm" 2 Hz" 1 s" 1 ). Here 1 nJy corresponds 
to m AB = 31.4. When needed, we have used the following 
relations to transform from Vega based magnitudes to AB 
magnitudes: K' = K' AB - 1.84, J = 3 AB - 0.82, I = l AB - 
0.42 and R = R AB - 0.17. 



2.2. Star formation rates 

A problem when using UV-luminosity densities to calcu- 
late the SFR is that these can be matched to almost any 
SFR, ranging from a strongly peaked to a flat, or even 
increasing, SFR at z > 1, by adjusting the assumed ex- 
tinction due to dust. By simultaneously using luminosity 
densities observed in different wavebands it is possible to 
break this degeneracy. Madau (1998) used three bands 
(UV, optical and NIR) in order to derive a SFR that can 
simultaneously reproduce the evolution of the different 
luminosity densities. He found a best fit for a universal 
extinction _Eb_v=0.1 with SMC-type dust. Despite cor- 
rections for absorption, the SFR still shows a pronounced 
peak at 1 < z < 2. This SFR is shown in Fig. [|. A peaked 
SFR is predicted in scenarios where galaxies form hierar- 
chically (e.g., Cole et al. 1994). Although commonly used, 
this SFR is not universally accepted and in Sect. 4 we 
discuss alternative SFR's. 

Core collapse SNe, which originate from short lived 
massive stars (ages <5 x 10 7 yrs), have an evolution 
that closely follows the shape of the SFR. Assuming an 
immediate conversion of these stars to SNe renders a 
multiplicative factor, k — SNR/SFR, between the SFR 
(Meyr^Mpc" 3 ) and the SNR (yr^Mpc" 3 ), where k de- 
pends on the IMF and the mass range of the progenitors. 
Here, we take 



rl25M Q 

JO.lMm 



M$(M)dM 



(1) 



Using a Salpeter IMF yields fc=0.0064. Note that, for con- 
sistency, the same IMF as assumed when deriving the SFR 
from the luminosity densities should to be used. Using a 
Scalo IMF, instead of a Salpeter IMF, decreases the con- 
version factor k between the SFR and the SNR by a factor 
2.6. A Scalo IMF, however, also increases the SFR as de- 
rived from the UV luminosities by a factor ~ 2, hence 
cancelling most of the effect. The reason is that the same 
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stars that produce the UV luminosities also explode as 
core collapse SNe. 

The lower mass limit for Type II progenitors is gen- 
erally believed to be 8-11 M© (e.g., Timmes et al. 1996). 
Here we choose 8 M©, a limit supported by e.g., Nomoto 
(1984). An increase to 11 M© would decrease k, and hence 
the SNR, by ~ 38%. Indications from especially the oxy- 
gen/iron ratio that the most heavy stars form black holes 
and do not result in SNe, justifies the use of 50 M© as 
an upper limit to the progenitor mass (Tsujimoto et al. 
1997). Other studies (e.g., Timmes et al. 1995), however, 
find that stars with masses down to ~ 30 M© may result 
in black holes. An upper limit of 30 M©, instead of 50 
Mq, decreases k by ~ 9%. 

The uncertainties concerning the origin of Type la SNe 
makes the relation between the rate of these SNe and the 
star formation more ambiguous. The predicted SNR de- 
pends on the nature of the SN progenitors. Yoshii et al. 
(1996) argue that the SNe la progenitor lifetime is prob- 
ably restricted to 0.5-3 Gyr. This range opens a possi- 
ble way to observationally distinguish between progenitor 
models. Ruiz-Lapuente & Canal (1998) find that the more 
short-lived double-degenerate progenitor systems and the 
long-lived cataclysmic-like systems should yield signifi- 
cantly different rates. 

In this paper we use as a standard model the SFRs de- 
rived by Madau (1998) to calculate the rates of both core 
collapse and Type la SNe. We use a universal extinction 
.Eb-v=0.1, together with SMC-type dust and a Salpeter 
IMF. In Sect. 4 we investigate how a higher extinction af- 
fects the counts, and if it is possible to use these counts 
to estimate the amount of dust. 

2.2.1. Core Collapse SNe 

The SNR used is shown in Fig. |l| This is the intrinsic rate 
of exploding SNe, irrespective of magnitudes and spectral 
distributions. The apparent magnitude of a SN at redshift 
z at a time t off its peak magnitude (i.e. t can be negative) 
in a host galaxy with inclination i observed in a filter / is 
given by 

m f (z,t,i) = M f (t)+n(z)+Kf(z,t)+Ag tf + < A u > .(2) 

Here, Mf(t) is the absolute magnitude of the SN in filter 
/ at time t relative to the peak of the light curve, (J,(z) is 
the distance modulus, 

fj,(z) = 51ogd L {z) - 5. (3) 
The luminosity distance, cIl, is given by, 

d L = «f^sinn{| n k |V> x 

x f [(l + zf(l + nMz)-z(2 + z)n A r 1/2 dz}. (4) 
Jo 

where Qj. = l — O^/ — Q\ and 'sinn' stands for sinh if Q k > 
and for sin if O fc < (Misner et al. 1973). If Q k = then 




Redshift 



Fig. 1. SNRs for the model with dust extinction 
Ss_y=0.1. The solid line represent the rate of core col- 
lapse SNe. This scales directly with the SFR, given by 
the right hand axis. Dashed, dotted and dash-dotted 
lines represent the rate of Type la SNe with times de- 
lays r=0.3, 1.0 and 3.0 Gyr, respectively. The rates of 
Type la SNe are normalized to fit the locally observed 
rates. 

the 'sinn' and the | fl k I 1 / 2 terms are set equal to one. 
For the standard CDM cosmology mainly used here, the 
distance modulus is given by fJ,(z)= 45.4 -51og(iJo/50 km 
s- 1 Mpc- 1 ) + 51og[(l+z) - (1 + zf/ 2 ]. Further, K f (z,t) 
gives the K-correction, A g j is the Galactic absorption, 
and < Aij > is the radially averaged absorption in the 
parent galaxy with inclination i. 

For Type Ib/c, plateau Type IIP and linear Type IIL 
SNe we use peak magnitudes given by Miller & Branch 
(1990). The magnitudes of the faint SN1987A-like SNe 
are not well known. Here we adopt the magnitudes given 
by Cappellaro et al. (1993), while magnitudes from Patat 
et al. (1994) are adopted for the Type Iln SNe. The mag- 
nitudes given by these authors are, however, not corrected 
for absorption. Adopting an average Eb-v =0.1 yields 
a mean As — 0.41 mag. Taking the effect of the albedo 
of the dust grains into account lowers the effective ab- 
sorption. We adopt an absorption As = 0.32 mag, which 
is consistent with the mean face-on absorption calculated 
by Hatano et al. (1998) in their models for Type II SN ex- 
tinction. We have here assumed that the absorption of the 
light from the parent galaxy and the SNe follow the same 
extinction law, implying that the SNe and the progenitor 
stars occur in the same environment, which to some extent 
may be incorrect. If the first core collapse SNe and their 
main sequence progenitors in a region sweeps away part 
of the shrouding dust, the absorption could be lower than 
calculated above. In Sect. 4 we return to this possibility. 
In Table [l] we list the corrected mean absolute B magni- 
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SN Type 


< M b ,q > 


a 


/ 


SNIb/c 


-18.44 


0.39 


0.23 


SNIIP 


-17.86 


1.39 


0.30 


SNIIL 


-18.38 


0.51 


0.30 


SN87A-like 


-15.43 


1.0 


0.15 


SNIIn 


-20.03 


0.6 


0.02 



Table 1. Adopted maximum absolute B magnitudes 
and dispersion in the AB-system. From Miller & Branch 
(1990), Cappellaro (1997), Cappellaro et al. (1993) and 
Patat et al. (1994). Values are corrected for extinction 
Eb-v=0-1, corresponding to Ab— 0.32. The intrinsic frac- 
tion of exploding SNc of the different types are given by 



tudes, as well as the adopted dispersion in this quantity 
for the different types at the peak of the light curve. 

For the evolution of the luminosity with time we use 
light curves from Filippenko (1997) for the Type Ib/c's, 
HP's, IIL's and the SN1987-like SNe. We have made some 
modifications to the light curve for the HP's to achieve 
a better fit to the observed, as well as theoretical, light 
curves presented in Eastman et al. (1994). For Type Iln's 
we use a light curve intermediate between IIP and IIL. 
This seems adequate in view of the data presented by 
Patat et al. (1994), who find that Type Iln SNe can have 
both linear and plateau shape, but there are also light 
curves in-between these two. 

The K-correction is calculated by assuming modified 
blackbodies for the spectral distributions of the SNe. At 
z^l even the I- band corresponds to the rest UV part of 
the spectrum. The exact spectral distribution in the blue 
and UV is therefore crucial. Unfortunately, this part of the 
spectrum is relatively unexplored even at low z. UV ob- 
servations of Type HP's are especially scarce. Only for the 
somewhat peculiar Type IIP SN1987A is there a good UV 
coverage (Pun et al. 1995). This showed already ~ 3 days 
after explosion a very strong UV deficit, similar to Type 
la SNe. This is a result of the strong line blanketing by 
lines from Fe II, Fe III, Ti II and other iron peak elements 
(e.g., Lucy 1987; Eastman & Kirshner 1989). Although the 
progenitors of typical Type HP's probably are red super- 
giants, rather than blue as for SN1987A, there is no reason 
why these ions should be less abundant than in SN1987A. 
On the contrary, because of the near absence of strong 
circumstellar interaction and similar temperature evolu- 
tion they are expected to have fairly similar UV spectra, 
as calculations by Eastman et al. (1994) also show. East- 
man et al. find that the UV blanketing sets in after ~ 20 
days when the effective temperature becomes less than ~ 
7000 K. This coincides with the beginning of the plateau 
phase. We therefore mimic the UV blanketing by a 2 mag- 
nitude drop between 4000 A and 3000 A, and a total cutoff 
short-wards of 3000 A, for the Type IIP and SN1987A-like 
SNe with T e //<i 7000 K. At higher T e ff we assume non- 



truncated blackbodies for these types. The time evolution 
of the spectra is modeled by changing the characteristic 
temperature of the blackbody curves to agree with the 
calculations by Eastman et al. At shock outbreak a short 
interval of energetic UV radiation occurs with 10 5 < T e ff 
^ 10 6 K. This only lasts a few hours, but may due to its 
high luminosity be observable to high z (Klein et al. 1979; 
Blinnikov et al. 1998; Chugai et al. 1999). It then cools to 
- 2.5xl0 4 K at - 1 day and further to <; 7000 K after ~ 
20 days. At the plateau phase the temperature is ~ 5000 
K. 

From the limited UV information of the Type IIL's 
SN 1979C, SN 1980K and SN 1985L (Cappellaro et al. 
1995), and the Type Iln SN 1998S (Kirshner et al., pri- 
vate communication), we model the spectra of the IIL's 
and Iln's by blackbody spectra without cutoffs. The most 
extensive coverage of the UV spectrum of a Type IIL ex- 
ists for SN 1979C (Panagia 1982). By fitting blackbody 
spectra to the optical photometry, we find that even after 
two months, when the color temperature is only ~ 6000 
K, there is no indication of UV blanketing. On the con- 
trary, there is throughout the evolution a fairly strong UV 
excess, consisting of continuum emission as well as lines, 
in particular Mg II 2800 A. The UV excess of these SNe 
is likely to be caused by the interaction of the SN and its 
circumstellar medium. This ionizes and heats the outer 
layers of the SN, decreasing the UV blanketing strongly. 
For the effective temperature as a function of time we use 
the fits for SN 1979C by Branch et al. (1981). The Type 
Ib/c light curve is taken from Filippenko (1997) and is 
assumed to have the same spectral characteristics as the 
Type la (described further below). 

The Galactic absorption, A g j, is taken to be zero in 
our modeling. This can easily be changed to other values 
of A g j. The term < Aij(z) > adds the absorption due 
to internal dust in the parent galaxy with inclination i, 
according to the adopted extinction laws. Gordon et al. 
(1997) show that the extinction curve for starburst galax- 
ies lacks the 2175 A bump, like an SMC-type extinction 
curve does, and shows a steep far-UV rise, intermediate 
between a Milky Way and an SMC-like extinction curve. 
The observed increase with redshift of the UV-luminosity 
originates mainly from starburst/irregular systems (e.g., 
Brinchmann et al. 1998). This implies that a major part 
of the core collapse SNe should be found in such galax- 
ies, and an SMC-type extinction curve should therefore 
be most appropriate when calculating the absorption of 
the SNe light. The difference between an SMC-type dust 
curve and a Milky Way-type extinction in the interest- 
ing wavelength range is small. This is especially true for 
SNe with a short wavelength cutoff in their spectral en- 
ergy distributions, but also a blackbody spectrum with 
T e / / ^ 7000 K drops fast enough at short wavelengths for 
the precise form of UV absorption in this region to be less 
important. 
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The dependence on inclination has been modeled by 
Hatano et al. (1998). The absorption closely follows a (cos 
behavior up to high inclinations. We approximate 
their results by < Ab,% > — 0.32[(cosi) _1 -l]. Also the ra- 
dial dependence of the absorption is discussed by Hatano 
et al. We simplify our calculations by adopting their ra- 
dially averaged value of the absorption. Dividing the host 
galaxies into different types increases the dispersion in ab- 
sorption. Finally, the extinction may depend on the intrin- 
sic luminosity and the metallicity of the host galaxy. This 
variation should to some extent already be accounted for 
in the observed dispersion of the peak magnitudes. 

The SNe are divided into the five different groups, with 
fractions, /, representing the intrinsic fraction of explod- 
ing core collapse SNe of the different types (i.e. irrespec- 
tive of magnitude). We estimate / by using the observed 
ratios of discovery, / f, s , given by Cappellaro et al. (1993), 
who find / o6s (IIP)~ / o6s (IIL), / o6s (Ib/c)~ 0.3/ o6s (II) and 
/ ohs (IIn)~ 0.2/ ofcs (II). 

The total number of Type II' s is the sum of Type IIP, 
IIL, 1987A-like and Iln. Adding the Type Ib/c yields the 
total number of core collapse SNe. Cappellaro et al. (1997) 
argue that the intrinsic fractions of Iln and 1987A-like 
should be /(Iln) = (0.02 - 0.05)/(II) and /(1987A-like) 
= (0.10 - 0.30)/(II). The IIP, IIL and lb Types have ap- 
proximately the same magnitudes, i.e. the ratio of dis- 
covery should be close to the intrinsic fractions between 
these. Combining these values and assumptions leads to 
our adopted intrinsic fractions in Table [l| We note here 
that, unfortunately, the rates of the different classes are 
affected by fairly large uncertainties. 

The observable number of SNe with different appar- 
ent magnitudes is calculated by integrating the SNR over 
redshift. The SNe are distributed between the different 
types according to their intrinsic fractions, and are placed 
in parent galaxies with inclinations between 0° to 90°. 
An important feature of our model is that the number 
of SNe exploding each year are distributed in time, and 
are given absolute magnitudes consistent with their light 
curves. With this procedure we obtain the simultaneously 
observable number of SNe, including both those close to 
peak and those at late epoch. In order to actually detect 
the SNe, at least one more observation has to be made af- 
ter an appropriate time has passed. In Sect. 8 we discuss 
the spacing in time between observations. 

For z > 5 we use a rate that is an extrapolation from 
lower z. This certainly simplifies the actual situation dras- 
tically. However, due to the large distance modulus and 
the decline in the rate at high z, the fraction of SNe with 
z > 5 is small, and hence the errors due to the uncertainty 
in the shape of the SFR at z > 5. Furthermore, SNe with 
a spectral cutoff at short wavelengths drop out at redshifts 
z ~ A e ///4000 -1, where A e // is the effective wavelength 
of the filter. In the R, I, J, K' and M' filters this occurs at 
z ~ 0.65, 1.0, 2.0, 4.5 and 10, respectively. This makes the 
contribution from SNe with higher redshifts insignificant. 



A caveat here is that lensing, as well as an early epoch of 
Pop III SNe, may cause a significant deviation from this 
extrapolation. These issues are discussed in Sect. 4 and 
Sect. 6. 

2.2.2. Type la SNe 

When calculating the number of Type la SNe we employ 
the same procedure as above. We use an extinction 
corrected peak magnitude Mb = -19.99 and dispersion a 
= 0.27, found by Miller & Branch (1990). The time delay 
between the formation of the progenitor star and explo- 
sion of the SNe is treated in a way similar to Madau et al. 
(1998a). Most likely, the progenitors are stars with mass 
3M Q < M < 8M (Nomoto et al. 1994). Stars forming 
at time t' reach the white dwarf phase at t' + AtMS, 
where AiMS=lO(M/M ) -2 - 5 Gyr is the time spent on 
the main sequence. After spending a time r in the white 
dwarf phase, a fraction -q of the progenitors explode as a 
result of binary accretion at t = t' + AtMS + t. The SNR 
at time t can then be written 

SNRQk) = v J* F SFR{t')dt'x 
x / 5(t-t' - M MS - r)<t>{M)dM, (5) 

J3M e 

where tp is the time corresponding to the redshift of the 
formation of the first stars, zp . Arguments from Yoshii et 
al. (1996) and Ruiz-Lapuente & Canal (1998) indicate 0.3 
<^ r <C, 3 Gyr. In the calculations we therefore use three 
different values for the time delay, r = 0.3, 1 and 3 Gyr. 
The r = 0.3 Gyr model approximately mimics the double 
degenerate case, while the r = 1 Gyr model resembles the 
cataclysmic progenitor model. With the additional r = 3 
Gyr the range is expanded to cover all likely models. 

The parameter r\ in Eq. (6) is introduced to give the 
fraction of stars in the interval 3 — 8 M Q that result in 
Type la SNe, and is determined by fitting the estimated 
SNR at z = to the locally observed Type la rates. For an 
alternative approach based on specific assumptions about 
the progenitors see J0rgensen et al. (1997). Results from 
local SN searches (Cappellaro et al. 1997; Tammann et 
al. 1994; Evans et al. 1989) give local rates of 0.12, 0.19, 
and 0.12 SNu, respectively (1 SNu= 1 SN per century per 
10 10 Lbq). Adopting a mean of 0.14 ± 0.06 SNu (Madau et 
al. 1998a), and using a local B band luminosity estimated 
by Ellis et al. (1996), leads to a local Type la rate of 1.3 
± 0.6 xl0~ 5 SNe yr _1 Mpc~ 3 . The normalization at z = 
yields an efficiency 0.04 < rj < 0.08, where the range 
is due to the fact that different values of r gives different 
r\ in order to reproduce the local rates. The uncertainty 
in the local SNR is equivalent to an additional spread 
in T) by a factor ~ 3. The normalization to the local rate 
leads to an uncertainty in the Type la SNR at all redshifts 
that corresponds to the uncertainty in the local rate. This 
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implies an uncertainty by a factor ~ 3 in the estimated 
Type la rates. 

The treatment of absorption in the case of Type la SNe 
is more complicated than in the case of core collapse SNe, 
and is therefore subject to larger uncertainties. The long 
time delay between the formation of the progenitors and 
the explosion unties the environmental link between these 
events. A (binary) star forming in a dusty starburst region 
may e.g. explode much later as a Type la SN in a dust-free 
elliptical galaxy. Observations show that Type la SNe do 
occur in both elliptical and spiral galaxies. However, the 
fact that the major part of local Type la SNe are detected 
in spirals, together with observations that indicate that 
the global fraction of ellipticals, or at least low-dust el- 
lipticals, seems to decrease at increasing redshifts (Driver 
et al. 1998) may justify our simplification of putting all 
the SNe la in spiral environments. Ignoring the elliptical 
parent galaxies leads to a slight overestimate of the ab- 
sorption, and a slight underestimate of the SN detection 
rate. 

Absorption in both the bulge and disk components 
of spiral galaxies have been calculated by Hatano et al. 
(1998). They find a somewhat smaller absorption for Type 
la SNe than for core collapse SNe. The Type la absorp- 
tion is also less dependent on the inclination of the parent 
galaxy. We use the disk component results of Hatano et 
al. in our model for the absorption of the Type la SNe. 
Fig. [I] show the intrinsic Type la SNRs for the three time 
delays used. Increasing the time delay shifts the peak of 
the la's towards lower redshift. The SFR (and the rate of 
core collapse SNe) peak at z ~ 1.55, while the Type la 
rates in the r = 0.3, 1 and 3 Gyr models peak at z ~ 1.35, 
z ~ 1.16, and z ~ 0.71, respectively. 

To describe the spectral energy distribution of the 
Type la SNe we use blackbody curves with a spectral cut- 
off at ~ 4000 A (e.g., Branch et al. 1983). The tempera- 
ture is set to 15 000 K around the peak, decreasing to 6000 
K after ~ 25 days (e.g., Schurmann 1983). We have com- 
pared the K-corrections in our model to those calculated 
from real spectra by Kim et al. (1996). The mean devia- 
tion in the R-band correction between our modified black 
body curves and the detailed calculations by Kim et al. 
is Am ~ 0.1 mag. This agreement justifies the use of the 
blackbody representation for the SN spectra. The average 
Type la light curve is taken from Riess et al. (1999). 

3. Results 

Using the model above we can now estimate the observed 
number of SNe per square arcmin down to some limiting 
magnitude, including corrections from extinction and the 
shift of the spectrum with redshift. 

As an illustration, we show in Fig. ^| the peak mag- 
nitudes (i.e. at ~ 10-15 days after explosion) for Types 
IIP, IIL, Iln and la as a function of redshift. We also show 
the magnitude of a Type IIP at the plateau phase (age ~ 



40 days). The magnitude for Type Ib/c's at peak follows 
the Type la's, except for an off-set of Am ~ 1.5 mag to- 
wards fainter magnitudes. The magnitude of SN1987A-like 
SNe at the peak resembles the Type IIP at the plateau, 
but with an off-set Am ~ 2.0 towards fainter magnitudes. 
Besides the standard flat Qm — 1 cosmology (hereafter 
SCDM), we also show the apparent magnitudes for an 
open cosmology (OCDM) with Dm = 0.3 and Qa = 0, 
and a flat, A-dominated cosmology (ACDM) with flm = 
0.3 and fl\ = 0.7. Note that the curves have a dispersion 
in magnitude according to values given in Table |. The 
dispersion is largest for the Type HP's; the peak mag- 
nitude of these may vary by more than one magnitude. 
Fig. H shows that the Type la's and the Type HP's at the 
plateau drop out of the I and K'-band at z ~ 1 and z ~ 
4.5, respectively, as a result of their UV-cutoffs. In con- 
trast, the UV-bright Type IIL's and Iln's stay relatively 
bright in the I band even at high z. 

3.1. Core collapse rates 

The solid lines in Fig. || show the number of predicted 
core collapse SNe per square arcmin in the R, I, K' and 
M' filters for different limiting AB-magnitudes. Because 
of the drop in the UV flux, bands bluer than R are of less 
interest for high redshifts. 

According to the specifications, NGST should have a 
detection limit of ~ 1 nJy in the J and K' bands and ~ 3 
nJy in the M' band for a 10 4 s exposure with S/N=10 and 
A/ A A = 3 (Stockman 1997). These fluxes correspond to 
AB magnitudes of Jab = 31.4, K' AB = 31.4 and M' AB = 
30.2, respectively. Using these limits, we find that the K' 
band yields the highest number of core collapse SNe. In a 
4x4 square arcmin field we predict ~ 45 simultaneously 
detectable core collapse SNe, with a mean redshift < z > 
= 1.9. 

In Fig. |we show the redshift distribution in all bands 
for mAB = 26 and thab = 31. For the lower magnitude 
limit one clearly sees the advantage of the infrared J, K' 
and M' bands when it comes to detect SNe with zj> 1. 
This is, of course, even more pronounced for thab = 31, 
with the M' band having the highest number of SNe at z> 
2. Note, however, that reaching Kab ~ 26 is considerably 
easier than Jab or K^ B ~ 26 from ground. The smaller 
detector size in the latter bands is also a severe limiting 
factor. 

For SNe with redshifts z > 5 we estimate ~ 1 SN per 
NGST field in the M' filter. The actual numbers of the 
high-z SNe are highly uncertain, since the SNR at z > 5 
is based on an extrapolation from lower redshifts. Also 
gravitational lensing may be important at these redshifts 
(see Sect. 6). However, independent of the actual numbers, 
we find that the M' filter yields a factor J> 2-5 higher 
counts compared to the K' filter at these redshifts. We 
also find that NGST should be well suited to detect SNe 
originating from a possible Pop III at redshifts J> 10. 
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Fig. 2. Magnitudes as a function of redshift for different SN types in the I-band (left panels) and the K'-band 
(right panels). Top panels show SCDM (Qm = 1, Qa = 0), middle panels show OCDM (Qm — 0.3, ft a = 0) and 
lower panels show ACDM (17m = 0.3, Ha = 0.7). Dotted line shows Type IIP at peak while the dash-dotted line 
shows the Type IIP at the plateau ~ 40 days after the explosion. Dashed, dash-triple-dot and solid lines show 
Type IIL, Type Iln and Type la, respectively. The magnitude for Type Ib/c's at peak follows the Type la's, 
except for an offset of Am ~ 1.5 mag towards fainter magnitudes. The magnitude of 1987A-like SNe resembles 
the Type IIP at the plateau, but with an offset Am ~ 2.0 towards fainter magnitudes. 



The R and I bands sample light with rest wavelengths 
short-ward of the peak in the blackbody curves at lower 
rcdshifts than the J, K' and M' bands do. Besides the 
drop in luminosity due to the spectral shape (an effect 



especially pronounced for SNe with strong UV blanket- 
ing), these wavelengths are affected by larger extinction. 
The large K-correction decreases the rates in the R and 
I bands, and few SNe with 1 are detected in I and 
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Fig. 3. Number of SNe per square arcmin that can be detected down to different limiting magnitudes in the 
R, I, K' and M' bands. Solid lines represent core collapse SNe. Dashed, dotted and dash-dotted lines represent 
Type la SNe with time delays r=0.3, 1 and 3 Gyr, respectively. Note that AB-magnitudes are used. 



R, even for limiting magnitudes ^ 29. Increases in these 
bands are instead caused by sampling SNe with fainter 
absolute magnitudes at z<, 1. 

In Fig. ^| we show the redshift distribution of core col- 
lapse and Type la SNe down to Iab = 27 and K' AB = 
31, respectively. The total number of core collapse SNe 
is shown as the solid lines in each panel. The lower solid 
line shows the fraction of these coming from Type IIL 
and Iln SNe, which lack UV-cutoff in their spectrum over 
the whole light curve. The Iab — 27 panel clearly demon- 
strates the dominance of these at high redshift. Besides 
these, only the fraction of Type HP's which are seen near 
the peak, before UV blanketing sets in, contribute to the 
SNe with 1.5. The K' band is far less sensitive to this 
effect, since the UV cutoff does not affect this filter at z<> 
4.5. When we compare the redshift distribution for K' = 
27 and K' = 31, we note that the mean redshift does not 
change much. However, the absolute number increases by 



a factor ~ 5. This is mainly caused by SNe below the light 
curve peak. Also the number of high redshift SNe with zj> 
2 increases by a large factor. 

Table || gives some examples of our estimates for the 
number of SNe for NGST, VLT/FORS and HST/WFPC2, 
all for an exposure of 10 4 s and S/N = 10. For VLT/ISAAC 
and HST/NICMOS the small field and the relatively 
bright limiting magnitudes, compared to NGST, do not 
result in more than ~ 0.01 SNe per field. These instru- 
ments are therefore of limited interest when it comes to 
detecting high-z SNe. 

So far we have discussed the number of SNe that are 
simultaneously observable during one search. To actually 
detect the SNe, additional observations are obviously re- 
quired. Preferentially, a series of observations of each field 
should be undertaken in order to obtain a good sampling 
of the light curves. This also leads to the detection of new 
SNe in the additional frames. The number of new SNe de- 



Tomas Dahlen & Claes Fransson: High-z Supernovae 



9 



—I 1 1 - 



—I 1 1 1— 



m AB =26 





—i — i — i — i — I — i — i — i — r 



2 J 4 5 6 
Redshift Redshift 

Fig. 4. Redshift distribution of core collapse SNe in the R, I, J, K' and M' filters for limiting magnitudes 
m/iB=26 (left panel) and m J 4s=31 (right panel). 



pends on the total length of the search, and the spacing in 
time between each observation. As an illustration, in an 
idealized situation where a field is covered continuously 
during one year, the detected number of SNe per square 
degree, with limit Iab = 27, is increased from ~ 260 in 
the first image, to ~ 1650 for the whole years coverage. 
This procedure is, of course, observationally unrealistic. 
Using 80 days between each observation (i.e. ~ five ob- 
servations for a years coverage) results in a total of 1150 
different SNe. Using 40 days instead (~ ten observations) 
gives ~ 1400 different SNe. Due to the time dilatation fac- 
tor (l+z), one observes relatively fewer new SN explosions 
at high redshifts compared to the first detection. 

With NGST, using limits as above, we estimate ~ 45 
SNe in each frame. The total number of different SNe in a 
field that is covered continuously during one year is ~ 68. 
This means that in addition to the ~ 45 SNe observed in 
the first field, only ~ 23 new SNe have exploded during 
the year. With three observations and 180 days between 
the observations, ~ 63 different SNe are detected. 

Compared to earlier estimates, our use of complete 
light curves during the whole evolution, as well as dis- 
tribution of the SNe over time, results in a larger number 
of SNe, as well as a realistic distribution over redshift for a 
given magnitude. For example, with the same SFR and ex- 
tinction, Madau et al. (1998a) predict ~ 7 SNe per NGST 
field per year in the range 2 < z < 4. Our calculations 
result in ~ 22. The difference is due to the fact that we 
use light curves covering the whole evolution, which allow 
us to include SNe at all epochs, instead of only those at 
peak, which is the case in Madau et al. 



3.1.1. Shock breakout supernovae 

Chugai et al. (1999) have noticed that the short peak in 
the light curve connected with the shock breakout may 
give rise to a transient event with a duration of a few 
hours. The possibility to observe this was pointed out al- 
ready by Klein et. al. (1979), although they concentrated 
mainly on the soft X-ray range. 

Based on a radiation-hydrodynamics code, similar to 
that of Eastman et al. (1994), Chugai et al. have calcu- 
lated monochromatic light curves for a Type IIP SN (or 
rather a scaled SN 1987A model) and a Type lib (specif- 
ically SN1993J). With a short time interval between ob- 
servations these SNe will be easily distinguishable from 
the Type la and Ib/c SNe, which have a rise time of t 
;> 20 days (at 1), and therefore only show a mod- 
est change in luminosity. A major problem in comparing 
their results to ours is that it is not discussed how they 
obtain their adopted intrinsic SNR's, although they ap- 
proximately agree with those used in this paper, as well as 
Madau (1998). They also neglect dust extinction in their 
calculations. 

Using this model, Chugai et al. find that two deep 
exposures, separated by ~ 10 days, result in 1.3 Type II 
SNe in the 6.8x6.8 square arcmin field of the VLT/FORS 
camera, using limit Iab — 28.2. In the calculations Chugai 
et al. assume that all SNe with z < 2 are detected with 
this limit. 

Using our hierarchical model, which gives approxi- 
mately the same SNR up to z ~ 2, and the same ob- 
servational set-up, as Chugai et al., our calculations result 
in ~ 0.27 SNe. The reason for our lower estimate is that 
Chugai et al. assume all Type II SNe to have the same 
steep initial rise as the Type IIP and lib. In our model we 




Fig. 5. Top panels: Observed SN distribution in redshift in the I filter, for limiting magnitudes 1,4.8=25 and 
Lab=27. Upper solid lines show the total number of core collapse SNe, while the lower solid lines show the part 
originating from Type IIL and Type Iln SNe, i.e the SNe without a UV cutoff in their spectra over the whole 
light curve. Dashed, dotted and dash-dotted lines show the number of Type la SNe for models with t=0.3, 1 
and 3 Gyr, respectively. Bottom panels: Same for K' filter with limiting magnitudes K^ B =27 and K^ B =31. 
Note the expanded redshift scale in the K' band, as well as the different scales of the y-axis. 



do not include any shock breakout for Types IIL and Iln, 
since the early time behavior of these SN types is not well 
known. If we do include all Type II SNe, our estimate in- 
creases to 0.66 SNe. The remaining discrepancy is mainly 
caused by the simplification Chugai et al. do by assum- 
ing that all SNe with z < 2 are detected, combined with 
the fact that they do not include dust extinction in their 
calculations. Nevertheless, this may be an interesting way 
of studying the shock breakout of SNe. Unfortunately, it 
may be difficult to estimate the temperature and lumi- 
nosity separately from this type of observations, because 
most of the observed evolution will be in the Rayleigh- 
Jeans part of the spectrum. Soft X-rays, as proposed by 
Klein et al., is here a better probe. 



3.2. Type la SNe 

Fig. U shows that in the R and I bands the number of core 
collapse SNe is comparable or larger than the number of 
Type la SNe for magnitudes fainter than 25 — 26. The 
exact crossing point depends on the life time of the SN la 
progenitors, as well as the Type la normalization at low z. 
In the K' and M' bands core collapse SNe tend to dominate 
at all magnitudes. The reason for this difference is that 
the Type la SNe have a higher effective temperature over 
a longer period than the Type II SNe. The optical to IR 
flux ratio is therefore higher for the Type la's. 

To illustrate the dependence on the progenitor life 
time, Fig. || gives the number of Type la SNe for r = 
0.3, 1 and 3 Gyr in the different filters. Fig. || shows that 
a change in the progenitor life time, r, introduces a non- 
negligible variation in the predicted number of observable 
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Table 2. Estimated number of SNe per field for different instruments. The limiting magnitudes are given for a 10 4 s 
exposure and S/N=10. 



SNe. For example, with an NGST detection limit, 
= 31.4, we predict ~ 8 SNe for the two low values of r, 
and ~ 5 for r = 3 Gyr. Observations in the I band with 
limits and field as for the VLT/FORS (see Table g) re- 
sults in 0.8, 1.1, and 1.8 SNe for increasing values of r. 
Counts may therefore seem like a useful probe to distin- 
guish between different progenitor models (Ruiz-Lapuente 
& Canal 1998). However, the uncertainty in the modeling, 
especially in the normalization of the Type la rates to the 
local value, makes the counts highly model dependent. In 
next section we show that this is further hampered by the 
additional dispersion introduced when considering alter- 
native star formation and extinction models. 

If both SN type and redshift information are available 
for the observed SNe, it may be possible to use the redshift 
distribution of the SNe to distinguish between progenitor 
scenarios. Figs. [j] and || show that the peak in the SNR 
moves to lower z as r increases. Also, the high redshift cut- 
off in the rates depends strongly on r. Fig. [j] shows that 
the rates decrease towards zero at redshifts z ~ 5.5, z ~ 
3.5 and z ~ 1.5, for models with r = 0.3 Gyr, r = 1 Gyr 
and r = 3 Gyr, respectively. To reach these redshifts, fil- 
ters unaffected by the UV cutoff must be used (i.e., A e //^ 
4000(l+z) A). Fig. | shows that the I filter is insensitive 
to Type la's at z J> 1.5, due to the spectral cutoff. Using 
the K' filter (lower panels of Fig. ||) makes it possible to 
sample all types of SNe up to z ~ 5, which is therefore 
most suitable for distinguishing different progenitor mod- 
els. As already mentioned, this requires a determination 
of both SN type and redshift. A discussion of methods and 
problems regarding this follows in Sect. 7. 

The only observational estimate of a Type la SNR at 
moderate redshift is by Pain et al. (1996). From a careful 
analysis, using realistic light curves and spectra, they find 
a Type la rate at z ~ 0.4 of 34.4^^ SNe yr" 1 deg~ 2 for 
21.5 < Rab < 22.5. Using the same magnitude interval, 
and counting the SNe exploding during one year, we find 
5-24 SNe yr" 1 deg~ 2 for r = 0.3 - 3 Gyr, where the 
lower number corresponds to the r = 0.3 Gyr model. The 
mean redshift of these SNe is in our calculation < z >^ 
0.35. These results seem to agree well, possibly favoring 
a high value of r. Note, however, that our estimated rate 



of Type la's is highly dependent on the normalization set 
by the local rate of these SNe (i.e., the efficiency parame- 
ter rj). We have already seen that the uncertainty in this 
normalization may be a factor ~ 3. 

3.2.1. Number of pre-maximum Type la SNe 

The number of simultaneously detectable SNe discussed 
above is a result of events over the whole light curve. Us- 
ing Type la's as standard candles for determination of 
flo requires observations at the peak of the light curve, i.e 
that a first detection is made at the rising part of the light 
curve. We estimate the number of such SNe by assuming 
that the comoving rise time is 15 days. Fig. || shows the 
number of Type la SNe down to different limiting magni- 
tudes before the peak of the light curve in the I filter for 
the three values of r. Also shown is the number of such 
SNe with z > 1 (the lower sets of curves) . An I band survey 
covering a one square degree field with limiting magnitude 
Iab ~ 27, will detect 16 — 31 pre-maximum Type la SNe 
(lower numbers for lower values of r), corresponding to ~ 
30% of the total number. About 5 — 8 of the pre-maximum 
Type la's have redshifts z > 1. These numbers are, again, 
sensitive to the local rate of Type la SN. 

4. Alternative star formation scenarios 

4-1. High extinction models 

There are several claims that an extinction £'s-y=0.1, 
as estimated by Madau (1998) and used here, is too low. 
Meurer et al. (1997) use the ratio between the far-IR and 
UV fluxes as a probe of the dust extinction in a galaxy, 
and find that the high redshift UV dropout galaxies may 
have a UV absorption of 2-3 mag. Using ISO observations 
of the HDF, Rowan- Robinson et al. (1997) estimate that 
only a third of the star formation is revealed by the UV- 
luminosity, with the rest shrouded by dust. Observations 
with SCUBA at 850 fim (Hughes et al. 1998; Smail et 
al. 1997) also indicate that a large fraction of the early 
star formation is hidden by dust. Several other estimates, 
based on observations of high- z galaxies (eg. Sawicki & Yee 
1998; Ellingson et al. 1996; Soifer et al. 1998), support a 
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Fig. 6. The number of simultaneously detectable Type 
la SNe on the rising part of the light curve for the I filter 
for three different values of t. Upper sets of curves show 
the total number of SNe, while the lower show SNe with 
z >1. 

higher extinction of Eb-v ~ 0.3. Observations of the far- 
infrared extragalactic background (Burigana et al. 1997) 
also seem consistent with a higher extinction. 

With an extinction of Eb-v ~ 0.3 it is likely that the 
observed peak in the SFR derived from UV-luminosities 
is illusive, and that the star formation history is compati- 
ble with a more constant rate, as in a monolithic collapse 
scenario (e.g., Larson 1974; Ortolani et al. 1995). In this 
model galaxy formation is thought to have occurred dur- 
ing a relatively short epoch at high redshift, zp > 5. This 
yields a SFR that increases from z = to z ^ 1.5 and then 
stays almost constant up to z J> 5. A SFR of this form may 
be compatible also with a hierarchical model, as shown in 
a recent paper by Somerville & Primack (1998). 

It should be noted that although there is much in favor 
for models with a flat SFR to zj> 5, they tend to over- 
predict the metallicity in the region z ~ 2 - 3, compared 
to estimates from damped Lyman-a systems (Blain et al. 
1999; Madau et al. 1998b; Pei et al. 1998). These mod- 
els also seem to over-predict the local K-band luminosity 
(Madau et al. 1998b). 

We have calculated the SFR and the resulting SNRs 
for this high dust scenario by adjusting the observed UV 
luminosities of Lilly et al. (1996) and Madau (1998) to an 
extinction Es_y=0.3. We also increase the local SFR in 
this model so that the observed B band luminosity den- 
sity of Ellis et al. (1996) is reproduced. The onset of star 
formation is set to zp = 5, while at even higher redshifts 
the SFR declines quickly. In Sect. 4.3 we comment on how 
higher values of zp affect the estimated rates. For the nor- 
malization of the type la's at z = 0, we assume that the 



increase in the intrinsic local B-band luminosity density 
should lead to an increase in the Type la rate by approxi- 
mately the same amount. In Fig. [?] we show this SFR and 
the intrinsic SNR for both core collapse SNe and Type la 
SNe, assuming r = 1 Gyr (dashed lines). 

4-2. Low extinction models 

The model above combines a high star formation at high 
redshifts with a high extinction in order to match the 
same observed luminosity densities. Another alternative 
scenario is an increased SFR at high z, but a low ex- 
tinction, as in the original hierarchical model. Evidence 
for such a scenario come from Pascarelle et al. (1998), 
who calculate the evolution of the UV luminosity den- 
sity up to z ~ 6 from ~ 1000 galaxies with photometric 
redshifts in the HDF. Taking the effects of cosmological 
surface brightness dimming into account and using a limit- 
ing surface brightness independent of redshift, they argue 
that the claimed UV luminosity density decrease at z > 2 
is mainly caused by a selection effect. Further, Hu et al. 
(1998) show that an increasing fraction of the total SFR 
at high redshift takes place in Lya emitters, and find that 
the fraction at z ~ 3 may be comparable to that derived 
from the dropout galaxies. Most of these objects would 
not show up in a survey such as the HDF, due to their 
low surface brightness. There is, however, no reason why 
SNe in these galaxies should not be detectable, especially 
since the authors argue that the extinction in these objects 
should be low. 

To mimic this scenario we have constructed a SFR 
with low extinction, Es_y=0.1, that increases like the hi- 
erarchical model up to z ~ 1.5, but stays flat at higher 
redshifts up to a formation redshift, zp — 5. The nor- 
malization of the Type la's at z = is the same as in the 
hierarchical model. This SFR and corresponding SNRs are 
shown in Fig. [j] as the dotted line. 

4--S. Results for alternative dust and star formation 
models 

We have repeated our calculations in Sect. 3 for the two 
additional scenarios described above. In the optical filters 
(probing rest-frame UV to optical), the increased SFR in 
the high dust model is mostly compensated for by higher 
extinction, resulting in observed rates that are only ~ 10% 
above those in the hierarchical scenario. The NIR filters 
(probing rest-frame optical to NIR) are less affected by 
extinction, which leads to a factor ~ 2 higher rates in the 
high dust model, compared to the hierarchical model. The 
estimates differ even more when comparing high redshift 
subsamples, i.e. when z > 2 is observed. This illustrates 
that the observed SNR may serve as an independent probe 
for the instantaneous SFR that is not subject to the same 
high uncertainty due to the unknown amount of dust ex- 
tinction as the UV-luminosity is. 
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Fig. 7. Intrinsic SNRs for the three different models. 
Solid lines show the hierarchical model with Eb-v — 
0.1, dashed lines show the high dust model with Eb-v 
= 0.3 and flat SFR at high z, and dotted lines show 
the model with Eb~v = 0.1 and a flat SFR at high z. 
Upper lines show the rates of core collapse SNe. These 
scale directly with the SFR, given on the right axis. 
Lower lines show the rate of Type la, for r = 1.0 Gyr 
model for each of the SFR models. 



For comparatively bright limiting magnitudes, that do 
not probe SNe above the proposed peak in the hierarchi- 
cal scenario, the two models with i?B_y=0.1 give, by con- 
struction, the same result. Observations must reach SNe 
at z ~ 2 before the rates start to differ. 

Except for the rare Type Iln's, a SN at peak magnitude 
at z ~ 2 has K' AB > 27, making NGST necessary for this 
type of observations. As an illustration, Table || shows 
to what extent rates of core collapse SNe are useful to 
constrain the SFR at high rcdshifts. For core collapse SNe 
with z ;> 2, NGST should detect ~ 3 and ~ 2 times higher 
rates for the high-dust-flat-SFR model and the low-dust- 
flat-SFR model, respectively, compared to the hierarchical 
model. At z> 4 these factors are ~ 5 and ~ 4. 

Due to the short time interval between formation and 
explosion in the case of core collapse SNe the choice of for- 
mation redshift enters only for observations that are deep 
enough to actually probe zp. Table ^ shows the variation 
of the estimated number of core collapse SNe with z > 5 
in the two models with flat SFR as Zf is changed from zp 
= 5 to zp = 7 and zp — 10. Note that more than 10% 
of the SNe in the M' band have a redshift > 5 in these 
models, compared to ~ 2% in the hierarchical model. 



Our calculations show a modest increase in the num- 
ber of Type la's for the high dust scenario, compared to 
the hierarchical scenario, at magnitudes brighter than ~ 
25. At magnitudes corresponding to z> 1, the differences 
are considerably larger. This is also true for the scenario 
with flat high-z SFR and low extinction, but the increase 
is modest, and it does not start before z ~ 1.5 As an 
illustration, Fig. || shows for the r = 1 Gyr model that 
for K' AB <, 31.4, the number of Type la SNe increases by a 
factor ~ 2 for the two models with flat SFR at high z, com- 
pared to the hierarchical model. Including the full range of 
t and different star formation scenarios, the NGST should 
detect between 5-25 Type la SNe per field in the K' band. 
For the ground based limit, Iab — 27, the number of Type 
la SNe increases by a factor 1.4-1.9 for the high dust sce- 
nario, and a factor 1.0-1.8 for the low dust scenario with 
flat SFR (smaller increase for lower values of r). 

If r is long, stars formed at early cosmological epochs 
may survive until low redshifts before ending as Type la 
SNe. Therefore, the increased star formation at high red- 
shifts in the two constant SFR models makes the rate of 
Type la SNe sensitive to the assumed formation redshift, 
zp. The choice of zp affects the SNR down to a redshift z, 
corresponding to a time t = tp + Atms + T - Increasing zp 
to > 5 increases the rates at zj> 2.7, zj> 2.0 and z^> 1.0 
for t = 0.3, 1 and 3 Gyr, respectively. For example, with 
r = 1 Gyr the high-z cutoff in the Type la rate occurs at 
z ~ 2.5 for the two models with flat SFR and zp = 5 (see 
Fig. 0). With zp = 7 and zp = 10 the cutoff in the rates 
moves to z ~ 3.0 and z ~ 3.4, respectively. This effect is 
larger for smaller values of r. 

As earlier mentioned, a further possibility that may 
increase the predicted counts of core collapse SNe is if 
the first SNe and their progenitors sweep away the dust, 
making the extinction lower for a large fraction of the 
SNe relative to the stars that produce the UV luminosity. 
We have studied this scenario by simply neglecting the 
extinction of the SNe in the two models with Eb-v =0.1 
(not included in Table ||) . The main difference occurs for 
the optical bands, where the effect of a lower extinction 
is largest. In the I band the counts increase by a factor ~ 
2, whereas the counts increase by a factor ~ 1.4 in the K' 
band for limits K Ag <^ 27. At fainter K' limits, i.e. reaching 
beyond the peak region, the differences decrease because 
a majority of all SNe is detected in both models, despite 
the different amount of absorption. For K' AB = 31.4 the 
increase is therefore only ~ 10%. 

5. Dependence on cosmology 

Besides the standard flat f2j\/ = 1 cosmology (SCDM) used 
so far, we have also studied two additional cosmologies; 
one open cosmology (OCDM) with SIm — 0.3 and f^A = 
0, and one flat, A-dominated cosmology (ACDM) with 17m 
= 0.3 and Oa = 0.7. The most obvious effect of a different 
cosmology comes from the luminosity distance, dh (Eq. 
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Fig. 8. Redshift distribution for core collapse and Type la SNe observed to K' AB =31 for the three different SFR 
models. Results for core collapse SNe are shown in the left panel, while the right panel shows the results for 
Type la SNe with r=1.0 Gyr. Note the different scales on the y-axis. 
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0.05 


0.05 


0.05 






K'=31 
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0.06 


0.02 


0.02 


0.02 


High dust, 


0.3 


M'=30 


6.6 


6.0 


3.6 


0.72 




0.75 


1.0 


flat SFR 




K'=31 


5.8 


5.2 
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0.29 




0.15 


0.20 


Low dust, 


0.1 
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3.2 


2.9 


1.9 


0.42 




0.44 


0.63 


flat SFR 




K'=31 


3.5 


3.1 


1.7 


0.25 




0.14 


0.20 



Table 3. Estimated number of core collapse SNe per sq. arcmin in the M' and K' band for different SFR and dust 
models. N(tot) gives the total number of SNe, while the other columns gives the number of SNe with a redshift above 
different specified values. The two models with flat SFR at high z has a formation redshift zp = 5, except in the two 
rightmost columns where the formation redshift is set to Zf — 7 and zp = 10, respectively. 



5). Changing from SCDM to OCDM and ACDM increases 
the distance modulus, making the high-z SNe fainter (Eqs. 
4 & 5), as seen in Fig. ||. The increased distance modu- 
lus in the OCDM and ACDM cosmologies also affects the 
intrinsic star formation rates used. A larger distance mod- 
ulus implies an increase in the absolute magnitudes of the 
galaxies from which the star formation rate, and hence 
also the supernova rates, are derived. As expected, these 
two effects, i.e fainter apparent SNe and an increased SNR, 
almost cancel when it comes to the observed number of 
core collapse SNe per square angle. 

The cosmology also affects the volume element, given 

by 

where cLm is the proper motion distance, <Im = (1+z) _1 c?l. 
The luminosity distance, cZl, is given by Eq. (5). The 
change in volume element affects the SNRs when ex- 
pressed in units number of SNe per Mpc 3 per yr. It does 



not, however, change the estimated rates of observed core 
collapse SNe expressed in number of SNe per square an- 
gle. This is due to the fact that the core collapse SNR is 
directly proportional the observed luminosity density of 
galaxies, which is calculated from the number of galaxies 
per redshift interval for a specific angle over the sky, an 
observational quantity independent of cosmology. 

When it comes to the Type la SNe, the volume ele- 
ment enters the calculations since the time between the 
formation of the progenitor star and the explosion of the 
SN may be a significant fraction of the Hubble time. This 
dependency on cosmology increases as the delay time t 
increases. The general trend is an increased SNR for the 
alternative cosmologies at high redshift. 

Fig. ^| show the estimated number per square arcmin 
of core collapse and Type la (using r = 1 Gyr) SNe in the 
I and K' filters down to different limiting magnitudes for 
the three cosmologies, using the hierarchical model for star 
formation. For the number of core collapse SNe the effect 
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Fig. 9. Number of SNe per square arcmin that can be detected down to different limiting magnitudes in the 
I and K' bands for three different cosmologies; SCDM (ft M = 1, O a = 0), OCDM (fi M = 0.3, Q A = 0) and 
ACDM (£Im = 0.3, Ha = 0.7). For the Type la rates we assume a delay time r = 1 Gyr. 



of changing cosmology is small for the reasons discussed 
above. Only at the very faintest magnitudes is there sig- 
nificant deviation between the models. At K' = 31, the 
estimated rates increase by a factor ~ 1.5 when changing 
from SCDM to ACDM. 

The estimated rates of Type la SNe differ at faint mag- 
nitudes by up to a factor two between the cosmologies. For 
t = 1 Gyr the intrinsic rates of the OCDM and ACDM 
models are higher than the SCDM at redshifts zj> 1. It is, 
however, necessary to reach faint magnitudes (K' ~ 28) 
to observe this increase in the total rates. Note, however, 
that a Type la at peak magnitude has m ~ 25 at z ~ 1. 
This means that using SNe seen at peak allows probing of 
the region where the rates start to differ between the cos- 
mological models at more moderate limiting magnitudes. 
Ruiz-Lapuente & Canal (1998) discuss the use of Type 
la's to distinguish between different cosmologies. In Sect. 
8 we comment on their results. 

The increased Type la rates at high z for the alterna- 
tive cosmologies means that the redshift cutoff moves to 
higher redshifts. This cutoff is also highly dependent on 
the delay time r (see Fig. [I]). It is, however, somewhat less 
dependent on the SFR. Therefore, if r is known, it should 
be possible to constrain the cosmology, even if the SFR 
is not accurately known. As an illustration, including all 
SFR scenarios, and using r = 1 Gyr (t = 3 Gyr), results 
in a cutoff in the Type la rates at z 2.5 - 3.5 (1.2 - 1.5) for 
the SCDM cosmology, while the drop is at z 2.9- 4.9 (1.6 - 
2.2) and z 3.2 - 5.3 (1.9 - 2.6) for the OCDM and ACDM 
cosmologies, respectively. The higher values in these red- 
shift ranges correspond to the hierarchical star formation 
scenario. 



6. Gravitational lensing 

In a recent paper Marri & Ferrara (1998) study the effects 
of gravitational lensing (GL) on high-z objects, in particu- 
lar Type II SNe, for a hierarchical model of galaxy forma- 
tion. For the three flat cosmologies they study (SCDM 
with Q M = 1, LCDM with Q M = 0.4, Q A = 0.6 and 
CHDM with n M = 0.7, tt u = 0.3), they find that there is 
at least a 10% chance that objects with z > 4 are magni- 
fied by a factor J> 3. To estimate the effects of GL on our 
results we have therefore used their magnification proba- 
bilities on our model (i.e. their Figs. 4 and 5). 

For the SCDM model, which yields the highest mag- 
nification, we find that the increase in the total number 
of core collapse SNe is only a few percent when using the 
NGST limits (M' AB = 30.2 and K' AB = 31.4). The effects 
are even smaller for the shorter wavelength bands. The 
reason for the small effect is that with the faint limits of 
the NGST almost all SNe up to z ~ 4 are detected even 
without magnification, and that the number of SNe with 
even higher z, where the magnification has largest effect, 
is relatively small. 

More interesting, for SNe with z >4 we find, when us- 
ing the NGST M' limit, an increase by ~ 20% of the pre- 
dicted counts. For SNe with z|> 9 the increase is ~ 40%. 
It should, however, be noted that estimates presented by 
Porciani & Madau (1998) give a much lower probability 
for substantial magnification compared to the Marri & 
Ferrara results used in the estimates above. The reason 
for this seems to be that Marri & Ferrara assume point- 
like lenses, whereas Porciani & Madau use a more realis- 
tic mass distribution characterized by singular isothermal 
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spheres. Therefore, the effects of gravitational lensing pre- 
sented here could be overestimated. 

7. Redshift determinations 

A major problem when comparing predictions from differ- 
ent scenarios to observations is the determination of the 
SN redshift. A direct spectroscopic determination with a 
resolution of > 100 is only possible for SNe more than 
two magnitudes brighter that the limiting magnitude, i.e. 
Iab ~ 25 for 8 — 10 m class telescopes and K' AB ~ 28.5 for 
the NGST. To reach fainter magnitudes the main alterna- 
tive is photometric redshifts of either the host galaxy or 
the SN. Photometric redshifts for galaxies have been dis- 
cussed extensively by e.g., Fernandez-Soto et al. (1999), 
Yee et al. (1998), Gwyn (1995) and Connolly et al. (1995). 
The fact that a large fraction of the star formation up to 
z> 1 occurs in dwarf galaxies, as well as the cosmologi- 
cal dimming, can make such a determination difficult. An 
alternative is to estimate a photometric redshift directly 
from the SN. A problem here is that the SN spectrum 
changes with both type and epoch. 

To examine this possibility we have determined broad 
band colors for the different types of SNe as function of 
epoch and redshift. As an example we show in Fig. |To| the 
color indices for a Type IIP SN as function of redshift at 
different phases. The spectra are taken from the synthetic 
spectra calculated by Eastman et al. (1994). These spectra 
assume LTE, but may nevertheless give a good impression 
of the qualitative evolution. 

At early time, before ~ 20 days the spectrum is a fairly 
smooth blackbody without a strong UV cutoff. The color 
indices do therefore not change dramatically with redshift. 
At later stages in the plateau phase the spectrum does not 
change much. An important aspect is that the UV cutoff 
at ~ 4000 A has now developed to its full extent, and the 
UV is essentially black. This is probably the most useful 
feature for identifying high redshift SNe photometrically. 
The extent of this UV drop may, as we discuss below, 
depend on the metallicity. The UV cutoff has a very pro- 
nounced effect on the optical color indices at z ^ 1, with 
strong increases in the B-I, V-I and R-I indices at succes- 
sively larger z. For z |> 1 the J-I and K-I, and finally K-J, 
are most useful due to the essential disappearance of the 
SNe in the optical. 

A problem with photometric redshifts of SNe, com- 
pared to galaxies, is that the colors, as we have seen, de- 
pend sensitively on the epoch. In addition, they depend 
on SN type. E.g., Type IIL's have less UV blanketing, 
while the Type la's have a rapid development of the UV 
cutoff. To break this strong degeneracy it is essential to 
have information about both SN type and epoch. It is 
therefore necessary to obtain reasonable light curves, i.e. 
a fairly large number, J> 5-10, observations of the field. A 
complete analysis can therefore be quite costly in terms 
of observing time. 



An alternative redshift method may be to use reason- 
ably well sampled light curves in combination with the 
cosmic time dilation. For Type la SNe one can safely as- 
sume a standard light curve. Although the absolute lumi- 
nosity can vary by a large factor, Type HP's have a fairly 
well defined duration of the plateau phase, which lasts ~ 
100 days. Also Type IIL's and Ib's and Ic's have reason- 
ably standardized light curves. From the observed light 
curve one can then, at least for SNe with z~>l, get an ap- 
proximate redshift within ~ 25% from a comparison with 
the low z light curve templates. However, since the light 
curve should be followed over a decline of about ~ 2 mag- 
nitudes in order to achieve a type specific light curve, the 
gain in depth by using photometry instead of multi object 
spectroscopy is marginal. The photometric accuracy also 
decreases for these levels and the actual limit may be even 
higher than two magnitudes. Therefore, in practice little 
is probably gained compared to direct spectroscopy. 

8. Discussion 

8.1. Uncertainties in the models 

A major source of uncertainty is the treatment of dust ex- 
tinction. In Sect. 4 we showed how different assumptions 
about the dust extinction affect the estimated rates. An 
underlying assumption in each of the calculations is that 
the same amount of dust affects the UV-luminosity from 
the high-z galaxies (used to calculate the SNR) and the 
light from the SNe. If the UV-luminosity originates from 
regions with dust properties that differ considerably from 
the regions where the SNe originates, an extra dispersion 
in the estimated rates should be expected. It is, however, 
difficult to estimate the uncertainties in the rates intro- 
duced by this, since the distribution of the dust within 
the galaxies at high z is poorly known. 

Other sources of uncertainties are the range of progen- 
itor masses in Eq. (2) and the choice of IMF (Sect. 2.2), as 
well as the dependence on cosmology (Sect. 5). Also, the 
distribution of the different types among the core collapse 
SNe influences the estimates. Changes in the fractions of 
the faint SN1987A-like or bright Type Iln's are most im- 
portant. Due to their low luminosity, the SN1987A-like 
SNe will be too faint to be detected by ground based tele- 
scopes for 0.5 (U4bS> 27 at peak). Doubling the fraction 
of these SNe from 15 % to 30 %, or, at the other extreme, 
setting the fraction to zero, changes the total number of 
detected SNe by ±18% at these redshifts. Changing the 
fraction of Type Iln's has a slightly different effect. By 
increasing /(Iln) from 2% to 4%, which is the upper limit 
proposed by Cappellaro et al. (1997), the total number of 
detected core collapse SNe at Iab ~ 24 increases by ~ 
17%. At these magnitudes it is mainly low redshift SNe in 
the steep part of the SFR curve that contribute, explaining 
the relative large effect. At fainter limiting magnitudes, an 
increasing fraction of normal luminosity SNe are detected 
at, and even beyond, the peak. The relative increase due 
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to the Type Iln's is therefore marginal. For Iab = 27 the 
number of SNe increases by ~ 4% and for K' AB ~ 31.4 the 
number increases by ~ 1%. 

Adding the uncertainties, we find that the counts of 
core collapse SNe may vary by a factor more than two 
due to insufficiently known model parameters. 

The estimated rates of Type la SNe are subjected to 
even larger uncertainties. Besides the factors which also 
affect the core collapse SNe, the rates of Type la's depend 
strongly on the assumed progenitor scenario, and are more 
dependent on cosmology. Also, the normalization at z = 
introduces an additional uncertainty by a factor ~ 3. 
Considering this, it seems unlikely that counts of Type 
la's can be used to set any constraints on the model pa- 
rameters. With additional information about the redshift 
distribution of the SNe it should, however, be possible to 
constrain either the progenitor life time or the cosmology. 
One of these parameters must, however, be independently 
determined, since there is a degeneracy between the both 
when it comes to the estimated redshift distribution of the 
SNe. To decrease the uncertainties involved, well defined 
searches of SNe at low redshift are highly desirable. 

8.2. Dependence on metallicity 

With increasing redshift the mean metallicity decreases, 
although the dispersion may be higher that at present. It 
is therefore interesting to discuss the consequences of a 
lower general metallicity. 

For Type la SNe, Kobayashi et al. (1998) argue that 
for one of the most likely progenitor scenarios, based on 
super-soft X-ray binaries, the necessary condition for a 
Chandrasckhar mass explosion may not occur for a metal- 



licity of Z <, 0.1 Zq, consistent with the decrease in the 
galactic [Fe/O] ratio at this metallicity. The physical rea- 
son can be traced to the peak in the opacity curve at ~ 
10 5 K, produced by iron. In the scenario by Kobayashi et 
al., this corresponds to a drop in the Type la rate at z ~ 
1.2. At higher rcdshifts Type la explosions are inhibited, 
since [Fe/H] <J - 1 here. Including the possibility of a dis- 
persion in the evolution of the metallicity leads Kobayashi 
et al. to conclude that a cutoff in Type la rates should oc- 
cur at z = 1 — 2. The Kobayashi et al. scenario should 
observationally be similar to the r = 3 Gyr model, with a 
turn-on at z ~ 1.5. Because the life time in their model is 
r = 0.6 Gyr, the peak in the SNR should occur at 1, 
rather than at z ~ 0.7, as in the r = 3 Gyr model. 

For core collapse SNe a lower metallicity can have sev- 
eral effects. First, line blanketing in the UV may decrease 
somewhat for the Type HP's. Since this is mainly an ion- 
ization effect, rather than an abundance effect, it is, how- 
ever, likely that this effect is small. This is partly con- 
firmed by the calculations by Eastman et al. (1994), who 
find only a marginal decrease of the UV blanketing as the 
metallicity is decreased from solar to a tenth of solar. For 
extremely low metallicity like for Pop. Ill stars with Z <> 
10~ 2 Zq, blanketing may, however, decrease significantly, 
although the Balmer jump will still be present. 

Secondly, the relative fraction between different SN 
types may change. In particular, the number of blue, com- 
pact supergiant progenitors similar to Sanduleak -69° 202 
may increase. This would then lead to a larger fraction of 
faint SN1987A like SNe, decreasing the number of observ- 
able core collapse SNe at high z. We emphasize that the 
exact reason for the blue progenitor of SN1987A is not 
fully understood (e.g., Woosley et al. 1999). 
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Finally, the mass loss process of the SN progenitor may 
depend on the metallicity. An example is the known de- 
crease of the mass loss rate with decreasing metallicity for 
radiatively accelerated winds (e.g., Kudritzki et al. 1987). 
In the red supergiant phase dust driven mass loss may be 
less efficient (Salasnich et al. 1999). The importance of bi- 
nary mass transfer may also depend on the metallicity. A 
change in the mass loss rate with metallicity would change 
the relative proportions between the different core collapse 
types. In particular, a decrease of the total mass lost is ex- 
pected to lead to a decrease in the number of Type IIL, lb 
and Ic SNe, while favoring Type HP's. In addition, a less 
dense circumstellar medium medium could then lead to a 
decrease in the ionization by the circumstellar interaction, 
and stronger line blanketing for the Type IIL and Iln's. 
Note, however, that the mass loss process even for local 
red supcrgiants in their final phase is poorly understood. 

8.3. SNe as probes of star formation 

With large ground-based telescopes, and especially with 
NGST, it should be possible to detect SNe up to high red- 
shifts, and to estimate the rates of both core collapse and 
Type la SNe. We have shown that because the rate of core 
collapse SNe follows the SFR, it should be possible to use 
observed rates of these SNe to constrain the SFR. As we 
have seen, a major problem is the influence of dust ex- 
tinction. In this respect we note that the NIR bands have 
the advantage of being less affected by dust extinction 
than the observed UV-luminosities. At high rcdshifts these 
bands correspond to the optical rest wavelength bands, 
and have therefore a factor of 2 — 3 lower extinction than 
the UV bands. Different star formation models may there- 
fore better be tested by using the K and M bands. Reach- 
ing high redshifts (z ^ 2), corresponding to K' AB ^ 27, 
increases the differences in the predicted counts between 
the star formation models significantly (see Sect. 4.3). 

The estimated difference in the redshift-integrated 
counts between the hierarchical star formation model and 
the high dust star formation model is a factor J> 2. This 
is about the same factor produced by the uncertainties in 
the modeling. It is therefore difficult to use these counts 
to probe star formation scenarios, unless the parameters 
involved are better known. What seems more feasible is 
to use rcdshift subsamples to constrain the shape of the 
SFR. As shown in Sect. 4.3, for models with a flat SFR at 
high z, we estimate a factor ^4 — 5 more SNe with z <; 
4. A major problem here is to determine the redshifts of 
the SNe (see discussion in Sect. 7). 

A further problem is if a large fraction of the star for- 
mation takes place in galaxies with very large extinction, 
like M 82 or Abell 220 with A v ~ 5 - 10 magnitudes 
(instead of Ay ~ 1, as for our high dust model). The 
difference between the optical and UV extinction is then 
of less importance, leading to a decrease in the estimated 
differences between the models. A similar large extinction 



may be indicated by the results of the ISO observations of 
some deep fields (e.g., Flores et. al. 1999). Far-IR observa- 
tions is then the most reliable way of deriving the true star 
formation rate. An alternative is to use some other source 
tightly coupled to star formation, but not affected by dust 
absorption. If gamma-ray bursts are related to some class 
of core collapse SNe (e.g., Type Ic's), they may be such 
class of objects (Cen 1998). Also radio observations may 
be interesting in this respect. 

8.4- SNe and nucleosynthesis 

The study of the nucleosynthesis by direct observation of 
SNe is naturally affected by the same problem as the star 
formation rate. Unless the dust extinction can be deter- 
mined reliably in an independent manner, the true num- 
ber of SNe is difficult to derive. In addition to this, the 
metallicity yields for SNe of different masses is non-trivial 
to derive even at low redshifts (e.g., Fransson & Kozma 
1999). Only for SN 1987A and a couple of other SNe has 
this become possible. The alternative to use theoretical 
yields from collapse calculations is obviously less satisfac- 
tory. The lower metallicity may also affect e.g. the mass 
loss processes, as discussed in Sect. 8.2. This may change 
both the progenitor structure and the upper and lower 
limits for the core collapse and Type la SNe, as well as 
the heavy element yields. In our view one of the most in- 
teresting goals for the observation of SNe at high redshift 
may be to observationally study the differences between 
the SNe in the early universe and those today. 

8.5. Spacing between observations 

An important aspect concerning the detection of SNe is 
the spacing in time between the observations. In order to 
detect the SN, the magnitude has to change appreciably. 
The interval is primarily dependent on the shape of the 
light curve. Near the peak, where the SN changes rela- 
tively fast, a comparatively short time is sufficient. This 
applies to searches where detection of SNe on the rising 
part of the light curve is the main objective (e.g., searches 
for Type la's for Hq). Core collapse SNe, which have a 
flatter decline of the light curve, need a longer spacing. 
This is especially true for the Type HP's, which in the 
plateau phase decline by <, 1 mag. Unless a SN can be 
detected (against the host galaxy) with this precision it 
will be missed. The limiting magnitude of the search also 
affects the necessary spacing. A deeper search results in a 
higher mean z of the observed SNe. Due to the cosmic time 
dilatation, the light curves of these SNe are stretched in 
time, implying that a longer interval between two observa- 
tions is needed, ~ 100(1 + z) days. Therefore, in order to 
detect these SNe a deep observation with the VLT requires 
an interval of > 100 days, while a corresponding NGST 
observation requires approximately a years interval. 



Tomas Dahlen & Claes Fransson: High-z Supernovae 



19 



8.6. Comparison to other works 

Apart from the studies by Pain et al. (1996), Madau et al. 
(1998a) and Chugai et al. (1999), which we have already 
commented on, there is a number of related investigations. 

Marri & Ferrara (1998) have studied of the effects of 
gravitational lensing of high redshift SNe. Using a Press- 
Schechter formalism and gravitational ray-tracing, they 
determine the magnification probability as function of red- 
shift for different cosmologies. We have already discussed 
the implications of their lensing results for our simula- 
tions in Sect. 6. Marri & Ferrara use these magnification 
probabilities to estimate the observed magnitudes at high 
redshift. The fact that there is a relatively large probabil- 
ity, <; 10% for a factor of three or larger magnification for 
z <: 4, means that even SNe as distant as z ~ 10 may be 
within the limits of NGST. When estimating the observed 
magnitudes they, however, assume that the light curve is 
described by a Type IIP light curve without any disper- 
sion in magnitude, although as we have seen, the Type 
HP's show a very large variation in luminosity. They also 
assume a fairly high temperature, ~ 25 000 K during the 
first 15 days, which is twice as high as the models by East- 
man et al. (1994) give. This is especially important for the 
high-z SNe, and, as Marri & Ferrara show, a lower temper- 
ature makes the detectability considerably more difficult. 
Marri & Ferrara do not attempt any discussion of expected 
rates of the high-z SNe. 

The effects of gravitational lensing is also investigated 
by Porciani & Madau (1998). They find, as earlier men- 
tioned, a considerably lower probability for a substantial 
magnification than Marri & Ferrara do. Porciani & Madau 
present I band counts for Type la and core collapse SNe, 
both including GL, and without lensing. These counts are 
presented as the number of SNe in different magnitude 
bins (21 < Iab < 27), seen at the peak of the light curve 
for an effective observation duration of one year. This leads 
to lower estimates for the observable number of SNe com- 
pared to our estimates, where we include SNe detected 
over the whole light curve. To compare our results we 
have calculated counts in the same units as used by Por- 
ciani & Madau. We find a fairly good agreement between 
the core collapse counts (deviation by a factor ~ 2), but a 
somewhat worse agreement between the Type la counts. 
It should be noted that expressing rates in units of an 
effective observation duration requires an idealized obser- 
vational procedure (as we have shown in the examples in 
Sect. 3.1). 

Ruiz-Lapeunte & Canal (1998) discuss the possibility 
of using R band counts of Type la SNe to distinguish dif- 
ferent progenitor scenarios. They find, similar to our es- 
timates, that models with long-lived progenitors result in 
higher counts than models with short-lived progenitors. 
To use this as a probe they note that it is necessary to 
know the SFR better than a factor 1.5. However, the un- 
certainty in the SFR seems, as we have shown, to be larger 



than this. It should therefore be difficult to use counts 
to determine progenitor scenarios. Additional information 
about the redshift distribution of the SNe is required. 

The same authors also estimate the effects on the 
counts for alternative cosmologies. They find that a flat 
A-dominated universe (ACDM) should result in higher 
counts of Type la SNe than a standard cold dark matter 
universe (SCDM). The difference between the cosmolo- 
gies start at m« ~ 24, and increases at fainter limiting 
magnitudes. A somewhat smaller increase in the counts 
is found for an open universe with zero cosmological con- 
stant (OCDM). 

Our results for different cosmologies agree with the 
general trend of Ruiz-Lapeunte & Canal. Using counts to 
distinguish between cosmologies, however, requires both 
that the SFR is well known, and that restrictions can be 
set on the progenitor life time. If this is not the case, 
the degeneracy between the different parameters involved 
makes a distinction between cosmologies very difficult. 

In an interesting paper Miralda-Escude & Rees (1997) 
discuss the possible detection of very high redshift core 
collapse SNe at z J> 5. By requiring that a metallicity ~ 
10~ 2 Z@ is produced at z ~ 5, they estimate a rate of about 
one core collapse SN per square arcmin per year above z ~ 
5. Our extrapolated hierarchical model gives a rate of ~ 
0.05 SN per square arcmin per year above z ~ 5. This may 
favor models with a flat SFR at high z, which result in <~ 
0.4 (0.7) SNe per square arcmin per year above z ~ 5 when 
using zp = 7 (zp = 10). However, the metallicity used by 
Miralda-Escude & Rees may be overestimated by an order 
of magnitude (Songaila 1997), leading to an overestimate 
of the SNR by the same amount. Also, the redshift before 
which the metallicity is assumed to have been produced 
affects the comparison. Using z ~ 3, instead of z ~ 5, 
decreases the estimated number SNe given by Miralda- 
Escude & Rees by ~ 30%. More important, integrating 
our rates for redshifts above z = 3, instead of z = 5 as 
done above, results in a number of SNe that is a factor <~ 
4 higher. Other uncertainties in the estimate by Miralda- 
Escude & Rees include the actual fraction of the baryonic 
matter which is enriched by the SNe. 

Miralda-Escude & Rees limit their discussion to Type 
IIP SNe, and do not attempt a detailed discussion of the 
observed rates. The observed magnitudes compare fairly 
well with our magnitudes in the K and M bands, but are 
brighter in the optical and near-IR bands. The main rea- 
sons for this is that they use a higher effective temperature 
and that they do not take into account any line blanketing 
in the UV, as our models do. As we discuss in next section, 
the low metallicity may decrease this effect. Apart from 
these caveats, the discussion by Miralda-Escude & Rees 
provides an important constraint at high redshifts. 

Gilliland et al. (1999) report on the discovery of two 
high redshift SNe in the HDF (see also Mannucci & Fer- 
rara 1999). One of the SNe has a probable host galaxy at 
z ~ 1.3 (photometrically determined) and is likely to be a 
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Type la, whereas the other SN has a probable host galaxy 
at z ~ 0.95 (spectroscopically determined) and is possibly 
a Type II. Gilliland et al. also make detailed estimates of 
the expected number of Type la and Type II SNe in a HDF 
like search. With limiting magnitude Iab ~ 27.7 they find 
that ~ 0.32 Type la SNe should be detected in a search 
consisting of the HDF together with an observation of the 
same area made two years after the HDF. Using the same 
cosmological model as Gilliland et al. (Om = 0.28, Oa = 
0.72, Ho = 63.3 km/s/Mpc 3 ) and r = 1 Gyr, we estimate 
~ 0.6 Type la SNe for a similar search. For a flat Q.m = 1 
cosmology we estimate 0.7 SNe. The main reason for the 
difference in results is that Gilliland et al. use a constant 
Type la SNR over the redshift range of interest (0 < z <> 
1.5), which is considerably lower than the mean value of 
our rates out to z ~ 1.5. 

For Type II SNe Gilliland et al. estimate - 1.2 SNe 
in HDF style search. This is in good agreement with our 
results, even though the modeling differs in many aspects. 
We estimate <~ 1.0 core collapse SNe for the cosmology 
used by Gilliland et al., and 1.3 SNe for a Q.m = 1 cos- 
mology. 

Considering the small statistics, both estimates are 
consistent with the discovery of two SNe in the HDF. 

Sadat et al. (1998) discuss the cosmic star formation 
rate, using a spectrophotometric model for different as- 
sumptions of the dust extinction. From this they calculate 
SNIa and core collapse rates, but do not translate these 
into directly observable rates. Their SFR is a factor J> 3 
higher than ours, which seems mostly to be due to the 
use of different factors when converting the observed lu- 
minosity densities to the SFRs. This also leads to higher 
SNRs (their Fig. 2). Sadat et al. also presents a case for 
Type la rates with a different normalization. These rates 
(their Fig. 3) agree better with our estimates at low red- 
shifts. At high redshifts the Type la rates differ more due 
to different modeling of these SNe. 

J0rgensen et al. (1997) attempt a calculation of the ab- 
solute rates of Type la, II and lb SNe from a population 
model. Although in principle appealing, this model de- 
pends on the uncertain scenarios for the progenitors of es- 
pecially the Type la's, as we have already discussed in this 
paper. Any estimates will therefore be sensitive to these 
assumptions. They also neglect the distinction between 
Type HP's and Type IIL's, which most likely originate 
from different progenitors. Further, J0rgensen et al. as- 
sume in the calculation of the observed magnitudes in the 
different bands as function of redshift, that the spectrum 
is characterized by that at the peak. As we have discussed, 
the spectrum and luminosity vary strongly with time. The 
most serious deficiency is in our view their neglect of the 
magnitude variation, as given by the light curve, which as 
we have seen, changes the observed rates by large factors. 
Their estimates of the observed rates are therefore highly 
questionable. 



9. Conclusions 

Observations of high redshift SNe are of interest for several 
reasons. First of all, one has through these a direct probe 
of the nucleosynthesis and star formation of the universe. 
In practice, there are several obstacles for a quantitative 
study of these issues. The fact that a large fraction of the 
star formation, and thus the SNe, may be hidden within 
optically thick dust can make it difficult to determine the 
total SFR and SNR accurately. This is certainly true for 
the optical bands, where we have found that the predicted 
total number of core collapse SNe with z < 1.5 is rather 
insensitive to the assumed star formation scenario, as long 
as the star formation is calculated to match the same ob- 
served luminosity densities, and the same extinction is as- 
sumed for the UV light from the galaxies and the SNe. 
Observations in the near-IR are less affected by this, and 
offers a clear advantage to the observations of the far-UV, 
as used for Lyman break objects. However, if a large frac- 
tion of the star formation occurs in highly obscured star 
burst galaxies, also the near-IR rates are severely affected. 
A further advantage of using SNe as star formation indi- 
cators is that they are insensitive to surface brightness 
selection effects. A complication when it comes to study- 
ing the nucleosynthesis is that the yields of the supernovae 
may vary with metallicity. 

An important motivation for searches of SNe at high 
redshift is that one can from this type of observations 
learn something about the SNe themselves when observed 
in a different environment. In particular, differences in the 
fractions of the various core collapse subclasses, their spec- 
tra and luminosities may give new insight into the physics 
of the SNe and their progenitors. 

The number of core collapse SNe that can be detected 
with NGST, with its expected limiting magnitude K' = 
31.4, should be ~ 45 per field in a 10 4 s exposure, assuming 
a hierarchical star formation scenario. The mean redshift 
of these SNe is < z > = 1.9. About one third of the SNe 
have z ~> 2. The high dust model results in total counts 
in the K' band that are a factor ~ 2 higher than in the 
hierarchical model. The estimated number of SNe with z > 
2 in the K' band for NGST is a factor ~ 3 higher than in 
the hierarchical model. The model with flat SFR at high 
z, but with low extinction, result in a factor ~ 2 higher 
number of SNe with z J> 2, compared to the hierarchical 
model, for the NGST limit. 

An important practical point is that in order to de- 
tect especially the Type IIP SNe at high z it is necessary 
to have a spacing between observations of ~ 100 days for 
ground based telescopes, and ~ 1 year for deep observa- 
tions with NGST. Shorter time intervals do not allow for 
the luminosity of the SNe to decrease by an amount nec- 
essary for detection. 

When it comes to the observed rates of Type la SNe, 
we find that these are highly sensitive to the star forma- 
tion modeling. This is due to the fact that the Type la's 
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are less linked to the environment where their progeni- 
tors were formed. The uncertainty in the life-time of the 
progenitors, combined with the sensitivity of the Type la 
rates to the onset of star formation in the models with a 
flat SFR at high z, contributes to the difficulty with using 
Type la counts as probes for either different star formation 
scenarios or progenitor models. This is further hampered 
by the fact that even the local rate is highly uncertain, 
and that this propagates to other redshifts through the 
normalization of the rates at z = 0. Therefore, accurate 
measurements of the Type la rates at low z are most de- 
sired. 

Precise measurements of the Type la rates at z J> 1 
could constrain the parameters to some extent. For exam- 
ple, in a given cosmology there is a high redshift cutoff in 
the la rates at an epoch that depends strongly on r, but 
is less dependent on the star formation scenario. In agree- 
ment with previous studies we find that counts of Type 
la SNe can be used as cosmological probes. This does, 
however, require that both the SFR and the unknown life 
time of the Type la progenitors can be determined inde- 
pendently. 

We predict the number of simultaneously detectable 
Type la SNe per NGST field to be - 5 - 25, depending on 
progenitor model and star formation scenario. Additional 
uncertainties widen this range even more. Of the simulta- 
neously observable Type la SNe, about 5% are on the ris- 
ing part of the light curve. For a ground based telescopes 
with limiting magnitude Iab ~ 27 we predict 75 — 400 
Type la's per square degree of which ~ 30% are on the 
rise of the light curve. 

A major technical aspect of our work is that we have 
tried to incorporate as much knowledge as possible about 
the theoretical and observational properties of the differ- 
ent classes of SNe. In particular, we have found that the 
spectral evolution is important for the magnitudes in the 
different bands. A striking example is the sensitivity to the 
UV cutoff, which for most SNe dominate the evolution in 
the optical bands. We have also found that a proper treat- 
ment of the light curve can change the predicted rates by 
factors of three or larger. An important source of uncer- 
tainty in these estimates are the local frequencies of SNe of 
different classes, as well as their distribution in luminosity. 
More extensive surveys with well defined selection criteria 
are therefore of highest priority for reliable predictions at 
high redshifts. 
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